Topological classification of vortex-core structures of spin-1 Bose-Einstein condensates 
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We classify vortex-core structures according to the topology of the order parameter space. We 
develop a method to characterize how the order parameter changes inside the vortex core. We 
apply it to spin-1 Bose-Einstein condensates and show that the vortex-core structures are classified 
by winding numbers that are locally defined in the core region. We also show that a vortex-core 
structure with a nontrivial winding number can be stabilized under a negative quadratic Zeeman 
effect. 

PACS numbers: 03.75.Lm, 02.40.Rc, 67.85.Fg, 67.30.he 



I. INTRODUCTION 

One of the salient features of superfluidity is the quan- 
tization of vortices. Due to the single- valuedness of the 
macroscopic wavefunction, the circulation of a scalar su- 
pcrfluid is quantized in units of h/M, where h is the 
Planck constant and M is the mass of the constituent 
particles. However, the situation changes in a very sig- 
nificant manner when the system has internal degrees of 
freedom. 

Since spinor Bose-Einstein condensates (BECs), 
namely, BECs with spin degrees of freedom, have been 
realized in ultracold atomic systems, a number of studies 
have been conducted to understand macroscopic quan- 
tum phenomena, including experimental 1 6] and theo- 
retical 043 studies on quantized vortices. In a spinor 
BEC, the circulation is quantized in units of a rational 
fraction of h/M, or not quantized at all. This is be- 
cause the spin and gauge degrees of freedom are coupled 
in a nontrivial manner. Examples include half-quantum 
vortices in the spin-1 polar phase 13 [ and 1/3-quantum 
vortices in the spin-2 cyclic phase 
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Vortices in spinor BECs have been extensivel y in ves- 
tigated and classified using homotopy theory jl7H2l| . 
However, the conventional homotopy theory does not tell 
much about the structure of a vortex core. For the case of 
a scalar BEC, the particle number density should vanish 
at the vortex core in order to avoid the phase singularity. 
We call such a vortex, i.e., a vortex accompanied by a 
density hole, a singular vortex. In contrast, for the case 
of a spinor BEC, the particle density does not have to 
vanish at the vortex core. Since the spinor BEC is de- 
scribed with a multi component order parameter, even 
when one component has a phase singularity, other com- 
ponents can fill the vortex core. As a whole, the particle 
density may become nonzero. We call this type of vortex 
a nonsingular vortex. A nonsingular vortex was realized 
in a spin-1 ferromagnetic BEC [g]. Recently, vortex-core 
structures have been numerically investigated for spin-2 
BECs [||- 

In this paper, we classify the vortex-core structure from 
the point of view of topology. Here, we define the vor- 



tex core and the core structure as follows. In general, a 
spontaneously broken-symmetry state is characterized by 
a coset space G/H, where G is the full symmetry of the 
system and H is a remaining symmetry of the broken- 
symmetry state (isotropy group). This coset space is 
called an order parameter manifold (OPM). We denote 
it as A4i ~ G/H in this article. In general, the classifi- 
cation of vortices is give n by calculating the fundamental 
group of the OPM |17H2l| . However, this classification 
is applicable only to a region which is far away from the 
vortex because, close to the vortex core, the kinetic en- 
ergy associated with the circulating current increases and 
the order parameter goes out of the OPM. We define 
the vortex core as the region in which the order param- 
eter leaves the OPM, and the term vortex-core structure 
is meant to represent how the order parameter changes 
in the vortex core region [23|. A topological method of 
classifying the vortex core structure was first proposed 
by Mermin et al. [24[ and Lyuksyutov (25[: this method 
tells us whether a vortex is singular or nonsingular. How- 
ever, there remains a question as to what the vortex-core 
structure is when the vortex is nonsingular. We address 
this question in the present paper by developing a general 
method to classify the vortex-core structure of a nonsin- 
gular vortex. Our method gives all topologically possible 
vortex-core structures. We then apply our method to 
spin-1 BECs and find that the vortex core can accommo- 
date not just one state but many different types of states 
which arc arranged in a concentric fashion. We char- 
acterize such a complex structure by introducing a lo- 
cal winding number along the concentric circle. Because 
the homotopy theory can enumerate possible states but 
cannot demonstrate their existence, we show that such 
vortex-core structures can indeed be realized by solving 
a timc-dcpcndcnt Gross-Pitaevskii equation numerically. 

The paper is organized as follows. In Sec. UH we re- 
view the classification of vortex singularities according 
to Refs. [Mill and apply it to spin-1 BECs. In Sec. Hill 
we generalize the method in order to classify the vortex- 
core structure. In Sec. IIV[ we apply our classification 
method to spin-1 BECs. We show that vortex cores are 
made up of concentric patterns in different states. In 
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Sec. El we verify that these vortex-core structures can 
indeed be realized as energetically stable configurations 
by numerical simulations under a magnetic field. Finally, 
in Sec. IVI1 we summarize our results. 



II. REVIEW OF THE CLASSIFICATION OF 
THE VORTEX CORE 

A mathematical scheme to classify a state at a vor- 
tex core was first proposed by Mermin et al. [24[ and 
Lyuksyutov [25| who applied the scheme to classify 
vortex-core structures in the superfluid Helium-3 and 
those in a liquid crystal, respectively. In this section, we 
briefly review their classification method and illustrate 
the method for the case of a spin-1 BEC. 

A. Classification of the vortex core 

We start with the energy functional given by 

E = J dr(e kin + e int ), (1) 

where ekm is the kinetic energy density and £i n t is the 
interaction energy density. Suppose that the interaction 
energy density can be decomposed into two terms, t\ and 
£2, as 

£int = ei + £2, (2) 

and that the energy scales of these terms are quite dif- 
ferent, say, | e x | *C The ground state in a uniform 
system is given by minimizing the interaction energy ej nt . 
Note that when the symmetry of the system is sponta- 
neously broken in the ground state, we cannot uniquely 
determine the ground state. Instead, we can construct a 
degenerate space of order parameters that minimize ei n t. 
For example, if the global U(l) gauge symmetry is spon- 
taneously broken, all states obtained by applying U(l) 
gauge transformations to a ground state also minimize 
the interaction energy. Such a degenerate space of order 
parameters is called an OPM, which is defined by 

M x := {.#| v .g G G/H 4 ,} ~ G/H^ (3) 

where tp is a ground state of ei nt . We define a full sym- 
metry group of system G, the action of which does not 
change the interaction energy, and an isotropy group of 
ip, which makes ip invariant under a group action: 

Jfy = {.9 G G\giP = !/>}. (4) 

We also introduce an enlarged OPM M2 as a space of 
the order parameters that minimize e 2 , where Mi and 
M 2 satisfy 

MiQM 2 - (5) 



In the presence of a vortex, the order parameter be- 
comes space-dependent, and the kinetic energy density 
is roughly given by ekm ~ h 2 n/(2Mr 2 ), where n is the 
particle-number density, M the mass of the particle, and 
r the distance from the vortex core. Hence, the kinetic 
energy density becomes negligible far from the vortex 
core and the order parameter in that region belongs to 
Mi- On the other hand, close to the vortex core, the 
kinetic energy density becomes comparable to the inter- 
action energy density. In the region of |ei| < e^m ^ | | , 
ekin dominates the weaker interaction energy ei and the 
order parameter belongs no longer to Mi but to an en- 
larged OPM M 2 . 

According to the homotopy theory [T3-EH , a topologi- 
cally stable vortex is labeled with an element of the fun- 
damental group iri(Mi). Since the OPM near the vortex 
core is enlarged to M 2 , the topological structure near the 
core is classified with 7Ti(A^2)- Here, we consider a loop 
I in the OPM Mi that cannot continuously shrink to a 
point in M 1 . The homotopy equivalent class of I charac- 
terizes a vortex. The loop / is embedded in the enlarged 
OPM M 2 . If the loop I shrinks to a point in M 2 , the core 
is nonsingular because the state in the core is specified by 
a point in M 2 . On the other hand, if the loop / cannot 
continuously shrink to a point in M 2l the correspond- 
ing vortex is topologically stable near the vortex core. 
Thus, the vortex can shrink to a point when we extend 
the enlarged OPM M 2 to the entire degrees of freedom 
of the order parameter, which includes a point at which 
the order parameter vanishes. Therefore, the vortex core 
is singular. The above statement can be expressed by us- 
ing an inclusion map from 7Ti ( M \ ) to ~k\ (M 2 ) . Because 
Mi is embedded in M 2 , there is an inclusion map from 
iti{Mi) to tti(M 2 ): 

m{Mi) ^ -k x {M 2 ). (6) 

For a given 7 G iti(Mi), if $(7) = l c , the vortex labeled 
7 is nonsingular, whereas if $(7) 7^ l c , the core becomes 
singular. Here, l c is an identity clement of tti(M 2 ). The 
singularity is characterized by an image of $. 



B. Application to spin-1 BECs 

We apply the above method to spin-1 BECs. The order 
parameter of a spin-1 BEC is given by 

V = (Vi^o^-ifeC 3 , (7) 

where ip m (m = 1,0,-1) describes a macroscopic wave 
function for atoms in the magnetic sublevel m. The 
mean-field energy functional for a uniform system is given 
by [HI3 

E= dr (e ki „ + eint) , (8) 
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where ekin and e; n t are respectively given by 

1 h 2 



(9a) 



(9b) 



Here, the first and second terms in Eq. (|9b[) represent 
the spin-independent and spin-dependent energy densi- 
ties, respectively, and n and \F\ are the local density 
and local magnetization, given by 



m— — 1 



F = - 

n 



\ 



mm"/" 



(10) 



(ii) 



respectively, with / = (f x ,fy,fz) being a vector of the 
spin-1 matrices given by 



(12) 




The coefficients of cq and c\ are given by 



Co 



Cl 



47rft 2 d 2 — Of) 



(13) 



M 3 ' M 3 

where as (S = 0, 2) is the s-wave scattering length for 
the spin channel S. In order for the system to be stable, 
the interaction coefficients should satisfy c > 0. 

The mean-field energy E is invariant under the U{1) 
gauge transformation and 5*0(3) spin rotation, namely, 
the full symmetry G of the system is given by 

G = U(l) 4> x SO(3) fl (14) 

where 4> and / stand for the gauge and the spin rotation 
symmetry, respectively. An element g £ G acts on the 
order parameter ■»/>, where g is represented by 

g = e^e- lS * a e- lS ^e-' lS ^ . (15) 

Here, </> describes the gauge degree of freedom, and a, (3, 
and 7 denote the Euler angles in the spin space. 

The ground state of the total energy is known to be 
a ferromagnetic (FM) state for c\ < and a polar state 
for ci > 0, where their representative order parameters 
are given by 



7/> F = ^(1,0,0) T (Cl <0), 



t/> P = A / -(1,0,1) ( Cl >0). 



(16) 
(17) 



These order parameters are invariant under the following 
operations: 

ff F = {(e^ ) e- i ^)| 7 e[0,27r)}, 

= ^(1W, (18) 
Hp = {(1, e~ lf yP), (e t7T ,e- lf ^e- lf yP)\P € [0, 2tt)}, 

= SO(2) fu x(Z 2 W, (19) 

Here, Hp represents the U(l) spin-gauge symmetry and 
Hp is constructed from a rotation symmetry around the 
y axis in the spin space and the spin-gauge coupled Z 2 
symmetry. The subscript <fi ± f z means the spin-gauge 
symmetry in terms of an operation (e 17 , e^^ 7 ). The el- 
ement of i?F(P) describes a group action on if) such that 

(e l0 ,e- l/ " 7 )^ := e ie e-^ 7 V V 0,7 e [0, 2tt), (20) 
where v = x, y, z. Therefore, the OPMs are given by [13, 

m 

G/Hp ~ (U(1U x SO(3) f )/U(l)t +f „ 

~ 50(3)0-/, , (21a) 
G/ffp * (C/(l) x SO(3) / )/(SO(2) / , x (Za)* +/ J 

(21b) 



([7(1)^ x S 2 )/(Z 2 ) 0+/ - . 



For real systems such as 8 Rb and 23 Na, the interac- 
tion coefficients satisfy cq 3> |ci|, which leads to 



c n 2 ^ jci|n 2 |J1|2 



2 2 
Thus, ei and e 2 are defined by 



ci"- 2 2 f _ c o" 2 

e l — — 7^1*1 ) £ 2 — — ^— ■ 



(22) 



(23) 



The enlarged OPM A^ 2 is given by minimizing e 2 , where 
the minimum is attained when the local density is con- 
stant. In this case, the enlarged OPM is given by 

M 2 =i S 5 . (24) 

To compare energy scales, we define characteristic 
lengths for eki n ~ ei, £2 as 



ft 



V2c raM ' 

h 

V 2 |ci|nM' 



(25) 
(26) 



where £o and £ s represent the healing length and the spin 
healing length, respectively. For cq 3> |ci|, we have 



Co < &• 



(27) 



The r dependence of the OPM for c 3> |ci| is reported 
in Table HI where r is the distance from the core. 
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TABLE I: The enlarged order parameter manifold in spin- 
1 BECs in the regions r < £o, £o iS r IS fs, and £s r - 
Here, r is the distance from the vortex core, and £o and 
£ s are defined by Eqs. ([5SJ and (J2HJ), respectively. In the 
rightmost four columns, we list the homotopy group for each 
enlarged OPM. Homotopy groups 7Ti, 7T2, and 7T3 character- 
ize vortices, monopoles (or two-dimensional Skyrmions), and 
three-dimensional Skyrmions in the three-dimensional space, 
respectively. The enlarged OPM for the superfluid helium-3 
is discussed in Ref [2(il |. 



Phase 


r 
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SO(3U- fz [27] 
x S 2 f ]/(Z 2 
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z 




z 


z 
z 



Since the fundamental group of the enlarged OPM is 
trivial, i.e., tti(S 5 ) = 0, elements of tti(A4i) are always 
mapped to an identity element: 

$: 7 ^l c , v 7 G7Ti(Mi), (28) 

that is, vortices are nonsingular. In fact, by the numeri- 
cal calculation, the energetically stable vortex-core struc- 
tures in spin-1 BECs have been achieved by mixing of the 
FM state and the polar state [l2j . 

III. GENERAL THEORY FOR CLASSIFYING 
VORTEX-CORE STRUCTURES 

As we have seen in the preceding section, whether or 
not a vortex is singular is determined by the inclusion 
map (|6]). Thus, it is natural to ask what the state inside 
the vortex core is when the vortex is nonsingular. To 
investigate the vortex-core structure, we generalize the 
classification method explained in the previous section. 

We consider an enlarged OPM which does not in- 
clude a zero point of the order parameter in order to 
characterize the nonsingular vortex. For such an en- 
larged OPM, the particle density is nonzero at any point, 
which implies that the enlarged OPM is homotopic to 
5 2s_1 = {ip\ J2n=i IVVi| 2 = const.} for a given order pa- 
rameter ip = (-01, • • • , ^g) € C s . Therefore, without loss 
of generality, we define the maximally enlarged OPM as 

V := S 2s -\ (29) 

The conventional OPM Mi, which is given by minimiz- 
ing the interaction energy Eq. ([3]) , is included in or equal 
to the maximally enlarged OPM V: 

MxCV. (30) 

In the case of spin-/ BECs, the maximally enlarged OPM 
is given by 

V = 5 4/+1 . (31) 



In the following discussion, we classify the core states 
and the structure of the order parameter inside the vortex 
core. We first decompose the maximally enlarged OPM 

V according to the symmetry property. Let G be the full 
symmetry group of the system under consideration, and 
we consider a G action on V. The G-orbit of ip 6 V is 
defined by 

M = {gi>\g € G}. (32) 

A set of all orbits in V is called an orbit space, which is 
denoted V in this paper: 

V:=V/G = m\^eV}. (33) 

In addition, we define an isotropy group of ip G V by 
Eq. Q. From the definition of the G orbit, elements of 
the G orbit [ip] are constructed by applying g € G to a 
representative order parameter ip. Since ip is invariant 
under the isotropy group H^,, we should subtract 
from G to obtain the G orbit. That is, the G orbit [ip] is 
equivalent to G/Hj,: 

m ^ G/H^. (34) 

When ip is the ground state that minimizes the interac- 
tion energy e; nt , G/H^, is identical to the OPM M.\. If 
ip and ip' belong to the same orbit, their isotropy groups 
are conjugate with each other, i.e., 

H^=gH v g-\ 3 g e G. (35) 

Note, however, that even when ip and ip' do not belong to 
the same orbit, their isotropy group may be conjugate to 
each other, which implies that different orbits may have 
the same property. Therefore, we categorize the G orbits 
according to the conjugacy class of their isotropy groups: 

(iP) :={[iP'] € V\H r =gH 4 ,g- 1 , 3 9 eG}. (36) 

Such a category is called a stratum. Then, the orbit space 

V is described as a union of (ip) : 

V=|JW- (37) 

ipev 

If V is compact, number of the strata is finite. More- 
over, since all isotropy groups for the order parameter in 
the same stratum are isomorphic to each other, can 
be replaced with a representative one Him . Then, from 
Eqs. (|3"3"|) and (|3"T|) . we obtain 

V ~ |J (G/H W x(iP)). (38) 

i>eV 

Next, we consider a map from a two-dimensional disk 
D 2 with radius R to V: 

O : D 2 -> V. (39) 

Here, we assume that the vortex is placed at the center 
of D 2 and the configuration of the order parameter on 
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vortex-core 
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Maximally enlarged OPM V 



Interaction energy 








0(R,9) GOPM G/H< VBS > 



FIG. 1: (Color online)We consider two-dimensional disk D 2 covering a cylindrically symmetric vortex configuration in a real 
system, where the core is placed at the center and the radius R is much greater than the healing length. Disk D 2 is mapped 
into the maximally enlarged OPM V, where the boundary dD 2 is mapped onto the conventional OPM denoted Q(R,6) with 
€ [0, 2n). With decreasing r from R to 0, the interaction energy increases because the loop enclosing the vortex gradually 
moves out of the OPM in order to shrink to a point. Finally, O(0, 9) is mapped onto the core state in the enlarged OPM. 



D 2 is axisymmetric. In this case, D 2 is parametrized by 
(r, 6>) with r G [0, R] and 9 G [0, 2tt). The boundary dD 2 
is defined by 

dD 2 = {(R,9)\9 G [0,2tt)}. (40) 

The map O is shown schematically in Fig. Q] Let 
"0BS G [ipBs] be a representative order parameter that 
minimizes the interaction energy ei nt , where BS stands 
for the boundary state; i.e., the OPM is given by 

M^G/H^y (41) 

When R is s ufficient ly large compared with the heal- 
ing length h/ y/2M\ti\, the image of the boundary dD 2 
should be a loop in the OPM, 

Q(R,6)eM^G/H^ Bs} , v 0G[O,27r). (42) 

The homotopy equivalent class of loop Q(R,9), which 
is an element of the fundamental group 7Ti(G/-H^ BS )), 
characterizes the vortex structure away from the core. 

To see the vortex-core structure, we need to investi- 
gate how loop 0(r, 9) moves in V as r goes to 0. We 
note that even when Q(R,9) is fixed in Mi, there are 
several possibilities for the map 0(r, 9) G V (0 < r < R), 
which is continuously deformed from Q(R,9). From the 
assumption that the configuration of the order parameter 
is axisymmetric about the vortex line, loop 0(r, 9) for a 
fixed r belongs to a G orbit [%p]: 

0(r, 9) G [V>] * G/H W , v G [0, 2tt), (43) 

Therefore, the loop <d(r,6) for every r is classified with 
the fundamental group iti(G/ Him). 

Since G/H^m is invariant as long as tp belongs to the 
same stratum, 0(ri,0) is nomotopic to 0(r 2 ,0) if [ip(r)] 
does not change the stratum for r*i < r < r 2 . In the 
case, the homotopy equivalence class of 0(fi, 9) and that 
of Q(r 2 ,9) belong to the same element of tti(G / H^). 
Accordingly, when r varies from R to 0, a nontrivial 



change of 0(r, 9) may occur when it crosses the bound- 
ary of strata. The nontrivial deformation of a loop can be 
characterized by a map between the fundamental groups 
7ri(G/iJ(^)) of the different strata. Hence, the vortex- 
core structure is characterized by a sequence of strata 
and by the map between their fundamental groups. 

In the following discussion, we define a map between 
the fundamental groups at the boundary of strata. We 
start by introducing the universal covering space of G, 
which we denote as G. The universal covering space is a 
simply connected space, namely, it satisfies 

G C G, (44a) 
tti(G)^0. (44b) 

If G is a simply connected space, we define G = G. As- 
sociated with the lift from G to G, the isotropy group 
H^) is also lifted to so as to satisfy 

G/H W ~G/H W . (45) 

Note that is not simply connected. When G is sim- 
ply connected, the fundamental group satisfies the rela- 
tions 

ni(G/H w ) £* MH W ) = H W /(H W ) , (46) 

where (H^ )o is a set of components that are connected 
to the identity element. Therefore, the map between 7Ti's 
is defined as a map between 's. We further note that 
the isotropy group of one of adjoining strata is bigger 
than or smaller than the isotropy group of the other, as 
proved by Michel [Tf|. It then follows that, when we 
consider a map from to , there are only two 

possibilities: (a) H/m C H/^n, and (b) Hlm D Hl^ix. If 
H( it ) C Hty>) (H w D H(i>')), the lifted isotropy group 
also satisfies Him C H/^>\ {Him 2 Hw))- Then, the 
map from Him to H/^i\ is defined as follows. 
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(a) For the case of H^ C H^i^ the map from H^ 
to Huj,i\ is defined by an inclusion map, 



i :H 



a H- a, 



v a 6 H 



(47a) 
(47b) 



(b) For the case of H^) 2 -^(i/>'}' we nrs ^ decompose 
Hup\ into a set of connected spaces {X{\ such that 



H, 



Xi n x 4 



f i,jeJ, (48) 



where J is a set of indices labeling Xi. When 
is a connected space, J includes only one element. 
We then define a map from each Xi to H^ as a 
constant map to an element b in the intersection of 
Xi and Hi, 



4° :Jsq->- fr w , &ejsqni%>, 

a i— > b, v a G A^. 



(49a) 
(49b) 



For each JQ, the number of constant ma ps i s equiv- 
alent to that of elements of Xi n H/^n |30l ] . 

The map between 7r (i?^) and tt (H^i) is induced by 
maps (j47aj) and (|49a|) : 



i*([a]) := [i(a)], y aeH w , 
(c^\([a]):=[c^(a)], v «a, 



(50a) 
(50b) 



where [a] denotes an element of /{H^)q. Note that 
map i is determined uniquely and its induced map is a ho- 
momorphism, whereas map c£ is not a homomorphism. 

The vortex-core structure is specified by giving all 
states between r = and r = R and assigning the wind- 
ing number to each state. Mathematically, this implies 
that the vortex-core structure is defined by the order of 
strata between r = and r = R and by the map between 
the lifted isotropy groups at each boundary between ad- 
jacent strata given in Eqs. (|50a|) and (|50bj) . Suppose that 
as r decreases, the order parameter moves on the strata 
along the following sequence: 



BS/ 



(A) {it) 



1-1/ 



(A)}- (51) 



Accordingly, the OPM changes from G/H^ Bg ^ to 
G/H^), and from G/H/^s to G/H/^ 2 \, and so on. 
Then, the vortex-core structure is defined as a sequence of 
maps between the fundamental groups of G/H/m, which 
is equivalent to the sequence 



toCHbs) 7T (-ffl 



MHi 



(52) 



where H^ := H/^ k \ (1 < k < I) is the lifted isotropy 
group of H/^ k } and J\. is a map defined in Eq. (|50a[) 



or 



(I50b[) . In order for loop 0(r, 9) to shrink in th< 
AG/Hi at r = 0, sequence (|52")) should satisfy 



_ shrink in the mani 

fold G/Hi at r = 0, sequence Q5 " 

Q;_l o ( _2 o • • • o Oi([/l]) 7^ L 
O; o o • • • o Oi([/l]) = [e], 



[e], 



(53a) 
(53b) 



where [h] is the element of ttq(Hbs) that characterizes 
the vortex under consideration, e is the identity element 
of Hi and the composition between maps is defined by 



n fc on fc _i([ft]) = fi fc (fi fc _i([/i])). 



(54) 



We call the final state of sequence, i.e. (^j), a core 
siaie and the sequence in Eq. (|52|) a core structure. In 
the following sections, we apply the above classification 
method to a spin-1 BEC. 



IV. 



VORTEX-CORE STRUCTURE IN A SPIN-1 
BEC 

A. Decomposition of 5 5 



From Sec. Ill B[ the maximally enlarged OPM is S 5 in 
a spin-1 BEC. In this section, we introduce the G orbit 
and the orbit space of S 5 . 

The order parameter in the spin-1 BEC is given by 
Eq. (J7J, which includes six independent variables. The 
full symmetry is given by Eq. (|14[) ; i.e., the energy func- 
tional is invariant under a gauge transformation and a 
spin rotation in the spin space. 

To derive the orbit space of S 5 , we consider a transfor- 
mation of variables. The physical quantities that charac- 
terize the state of the spin-1 BEC are the number density 
n and the magnetization \F\. It can be shown that these 
quantities, overall phase </>, and Euler angles a, j3, and 
7 completely specify the spin-1 order parameter. In gen- 
eral, the state of the spin-1 BEC is described by 



o 



2 



(55) 



where 0,a,7 € [0,2tt), @ € [0,tt], and < \F\ < 1. 
The derivation of Eq. (|55[) is shown in Appendix [X] 
Equation (|55j) tells us that the orbit space is spanned 
by n and \F\. When n is a constant, the total de- 
grees of freedom of the order parameter is equal to those 
of the five-dimensional sphere S" 5 , which corresponds to 
the maximally enlarged OPM. Thus, the orbit space is 
parametrized by the local magnetization < |.F| < 1, 
and, hence, given by 



V=[0,1]. 



(56) 



The orbit space V is decomposed into three strata: 

(ifr) - {0}, (57a) 
(V>c 2 )^(0,l), (57b) 
(#) - {1}, (57c) 
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G/Hp 



G/Ifo 



G/H F 



G-orbit 



[Oj 



(0,1) 



~° • orbit space 
{1} 



FIG. 2: Decomposition of S . By Eq. ([61}, S" is homeomor- 
phicto (G/Hp x {1}) U (G/Hc 2 x (0,1)) U(G/tf P x {0}). The 
orbit space is given by [0, 1], where [0, 1] represents the ampli- 
tude of local magnetization. The orbit space is decomposed 
into three strata: {0}, (0,1), and {1}. The corresponding G 
orbits are given by G/Hf at {0}, G/Hc 2 at (0, 1), and G/Hp 
at {1}. 



corresponding to the polar state (\F\ = 0), the Ci state 
(0 < \F\ < 1) [3l|, and the FM state (\F\ = 1). The 
order parameter and the isotropy group for the FM state 
are given by Eqs. (|T5|) and (|T5|) . respectively, and those 
for the polar state are given by Eqs. (fT7)) and (|T9|) . On 
the other hand, the C 2 state is defined partially 
polarized state whose representative order parameter is 
given by 



4>C<2 




(0<|F|<1). 

(58) 



The C 2 state is invariant under a it rotation about the 
z axis together with a tt gauge transformation, i.e., the 
isotropy group of the C 2 state is given by 



H C2 = {(1,1a), (e^e"*'-*)} S (Z 2 )^. 



(59) 



Here, 1„ denotes an n x n unit matrix. Note that the C 2 
state cannot be the ground state of the system in the ab- 
sence of a magnetic field (the magnetic field specifies not 
only the magnitude but also the direction of the sponta- 
neous magnetization [Hj]). In other words, contrary to 
the FM and polar states which become the ground states 
for ci < and c\ > 0, respectively, the C2 state involves 
excitations such as vortices. 

By the isotropy groups Hp and Hp , OPMs for the FM 
and the polar states are given in Eqs. (|21a[) and (I21bl) . 
respectively. By Eq. (|59|) . the coset space for the C2 state 
is given by 

G/H C2 s (C/(l) x SO(3) f )/(Z 2 )^ +fz , (60) 

for fixed 

From the above discussion, we find that S 5 can be 
decomposed in the form of Eq. (|38p. When the local 
density is constant, S is decomposed as follows: 

S 5 ~ (G/Hp x {0}) U (G/H Co x (0, 1)) U (G/H F x {1}), 

(61) 

where [0, 1] represents the amplitude of local magnetiza- 
tion |.F|. We show the decomposition of S 5 graphically 
in Fig. H 



B. Construction of maps 



We calculate the fundamental group of each state and 
derive the map i and 4° defined in Eqs. (f47b|) and (|49b)l . 
respectively. We begin by calculating elements of the 
fundamental group using Eq. (|46[) . Since the symmetry 
group G = U(l)<f, X 5*0(3) f in the spin-1 BEC is not sim- 
ply connected, we lift it to its universal covering space: 



G 



1+ x SU(2) f . 



(62) 



Associated with the lift, the isotropy group H q (q = 
F, P, C2) is also lifted to H q . As a result, Hp, Hp, and 
Hc 2 are given by 



Hp := {(x, ±e 
Hp := I (7m, e~ 



n G 



H, 



(•■> 



{< 



7m, ±e 



n € Z 



}■ 



(63a) 

Z, pe [0,4tt)|, 

(63b) 

(63c) 



where a u , (y = x,y,z) are the Pauli matrices. The 
derivations of Eqs. (|63a|) . (|63b[) . and (|63c|) are given in 
Appendix [Bj Here, we add the factor 7r in front of n € Z 
for the sake of convenience. In the following discussion, 
we describe each lifted isotropy group by the representa- 
tion of Eqs. (|o3a l> -([o3c ) . 

In the FM state, a set of components connected to 
(0, 1 2 ), (#f)o, is given by 



(Hp) = {(x,e-^ x )\xe 



(64) 



Hence, the zeroth homotopy group of the FM state is 
given by 

MHf) = Hp/(Hp) = {(0, l a ), (0, -1 2 )} 

= Z 2 . (65) 

The nontrivial element (0, — 12) represents a spin vortex. 
The structure of Hp is illustrated schematically in Fig. 
[31(a). 

For the polar state, (-ffp)o is given by 

(H P )o = {(0,e-^)|/?e [0,4tt)}, (66) 

from which the zeroth homotopy group is given by 

7t (#p) = Hp/(H P )o 

= {{2nir, 1 2 ), ((2n + 1)tt, eT l ^)\n e Z} 
^ Z. (67) 

Hence, vortices in the polar state are classified by an 
integer: (2-im, 1 2 ) describes a gauge vortex, whereas 
((2n + l)7r, e~ l ^ v ) represents a half-integer vortex (or 
Alice vortex). The structure of Hp is shown in Fig[3](b). 

For the C2 state, since Hq 2 is a discrete subgroup of 
G, (-Hm)o is composed of the identity element alone; 



( J ffc 2 )o = {(0,l 2 )}. 



(68) 
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(a) Hf (b)Hp (c)H C2 




FIG. 3: (Color online) Structure of isotropy groups (a) Hp, (b) Hp, and (b)Hc 2 in the universal covering group G. The 
vertical axis shows the transformation of SU(2)f represented by e~ l ~ p with /3 G [0, An), and the horizontal axis represents 
the transformation of R^, x SU(2)f given by (x, e~ 1 ^~ x ) with x G R. (a) Hp is composed of two connected spaces, X + = 
{(x, e~ l ~^ x )\x G R} and X- = {(x, —e~ 1 ~^ x )\x G R}, which correspond to the upper and lower horizontal lines, respectively. 

(b) Hp is composed of X„ p = {(nn p , e~ l ^~ 1Trlp e~ x ^~^)\l3 G [0,47r)}, where n p G Z. Each X Up corresponds to a vertical line. 

(c) Hc 2 consists of discrete elements X( + ,„ p ) = (ivn p , e~ l ^i~ nnp ) and -X"(_ jnp ) = (nn p , — e~'^~ nrlp ) with n p G Z, which are 
represented by filled circles. The thick (red) lines represent the nontrivial loop in the coset space G/H q (q — F,P,C2). Since 
all elements of H q are identified in the coset space, a loop in G/H q is described by a segment that connects the element of H q . 
From this figure, we can understand nontrivial loops such as (a) a unique loop, (b) an infinite number of loops labeled by single 
integers, and (c) an infinite number of loops labeled by two integers. 



Therefore, the zeroth homotopy group is isomorphic to 
Hc 2 - 



= {(27rn,±la),((2n+l)7r,±e- 



nezj. 
(69) 



Comparing Eq. (JSSJ) with Eqs. <JBS]) and flST]), ir (H C2 ) 
includes the elements of ttq(Hv) and those of iro(Hp). 
These structures are important to characterize the 
vortex-core structures as discussed in Sec. IIV CI The 
structure of Hq 2 is shown in Fig. [3] (c). 

Let us construct the map between 7To's of each lifted 
isotropy group. The lifted isotropy groups Hp, Hp, and 
Hc 2 satisfy 



Hc 2 C Hp, Hp. 



(70) 



Since ("0f(p)) does not adjoin (?/>p(f)), we have to pass 
through the C2 state. Therefore, we need to define only 
four types of maps: (a) from Hp to Hc 2 , (b) from Hq 2 
to Hp, (c) from Hp to Hq 2 , and (d) from Hq 2 to Hp: 



„(±) 



ip 



Hp , " > H C2 < K) > Hp. 
if Cn f 



(71) 



Each map is denoted c„, p ', ip, Cwt\ and ip, and defined 
as follows: 

(a) ch'J: Since Hp D Hq 2 and Hp consists of two 
connected spaces X + = {(x, e~ l ^~ x )\x G M} and 
X_ = {(x, — e~ z ~z~ x )\x G K}, the map is defined by 



Eq. (|49a| . i.e., 

c£J : Hp -> H C2 , (72) 

(x,±e- l ^ x ) h> {-Kn p ,±e- l ^ np ), (73) 

where x G M and n p G Z specifies an element of 
Hc 2 - 

(b) ip: Since Hq 2 C Hp, the map is given by the in- 
clusion map 



*f : H C2 —> Hp, 



(74) 



(iin p ,±e- l ^™ p ) H- (irn p ,±e- l ^ np ), n p G Z, 

(75) 

where ip is determined uniquely. 

(c) ch s p ^: Since Hp D Hc 2 and Hp consists of in- 
finite connected spaces labeled by an integer n p : 
X np = {(2Trn p ,e- l ^P e -^ 2 ™ p )\f3 G [0,4tt)}, the 
constant map is given by. 

c^:Hp^H C2 , (76) 

(27m p ,e- I ^<V l * 27m ») 1 — ^ (2irn p ,e-^ 2 ™f e- 1 ^ 2 ™"), 

(77) 

where (3 G [0,47r) and nf G {0, 1} specifies an ele- 
ment of Hc 2 - 

(d) ip: Since Hc 2 C Hp, the map is given by the in- 
clusion map 

ip : H C2 -> Hp, (78) 
(Tm p , ±e~ iS i-™p) 1 y (nrip, ±e~ lZ f™ p ), (79) 
where ip is determined uniquely. 
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(a) polar core vortex 
H F 



(0,1 2 ) 



(P,-W 

(b) ferromagnetic core vortex 



Hp 



H C2 



(0,l 2 ) 



(0,-l 2 ) 



(0,1 2 ) 



(0,-l 2 ) 



H C2 



Hp 



(0,1 2 ) 



(<=„«)* 



(ip)* 



(0,-l 2 ) 



Hp 



Hc 2 



(0,1 2 ) 



(0,-l 2 ) 



H,- 



(ip)* 



(0,1 2 ) 



(04 2 ) 



(0,1 2 ) 



► 



(i F )* 



(0,-l 2 ) 



(0,-l 2 ) 



FIG. 4: (Color online) Deformation processes for (a) the polar-core vortex with n v = and (b) the FM-core vortex with n p = 1 
in the FM state, (a) The vortex, which is topologically stable in the FM state, shrinks to a point in the polar state. In the 
intermediate region, we run through the C 2 state, (b) The vortex shrinks to a point in the FM state. In the intermediate 
region, we pass through the C 2 state and the polar state. The polar state has a nontrivial winding. 



Using these maps c Up , «f, chf , and ip, the maps be- 
tween 7To's are induced as follows: 

(eg)). : 7t (#f) -> MHc 2 ), (80a) 

(i F )* : MHc 2 ) ->■ 7r (#p), (80b) 

(cg^)). : 7T (flp) TTofe), (80c) 

(i P ), : ^o(#c 2 ) -> 7t (Hp). (80d) 

The definition of the induced map is given by Eqs. (|50a[) 
and (|50b| . 



C. Classification of a vortex-core structure 

1. FM state 

For the case of c\ < 0, the boundary state becomes the 
FM state, 

</>bs = V>f = v^(l,0,0) T (81) 

The nontrivial element of tto(Hf) is given by [(0,— 12)]. 
Here, [• • • ] represents the homotopy equivalent class. In 
what follows, we derive all possible vortex-core structures 
with the ferromagnetic boundary state. First, by using 



the map (c )*, we obtain 

(ip),o(4- ) ),([(0,-l 2 )]) = (ip),([(0,-l2)]) 

= [(o,-i 2 )] 

= [(0 ) la)]G7r (-ffp) ) (82) 

where (0, 12) is the identity element of H q (q = F, P, C2). 
For the last equality in Eq. (|82p . we use the homotopy 
equivalence relation in £/p, 

(0 J l 2 )~(0,-l 2 )G(Fp)o. (83) 

The deformation process in Eq. (|82p is illustrated 
schematically in Fig. 2] (a). Equation (|82l) shows that 
there exists a vortex whose core structure is given by 

7r (i? F ) MHc 2 ) MHp), (84) 

where sequence (|8"4"|) satisfies Eqs. (|53a[) and (|53b|) ; i.e., 

(4- ) ),([(0,-l 2 )]) = [(0,-l 2 )], (85a) 

( l p)*°( C ^ ) )*([(0,-l 2 )]) = [(0,l 2 )]. (85b) 

According to Eq. (fSTj) . we describe the order of strata of 
Eq. ((SI} as 

{(^f) -> <</>C 2 ) -> (Vp)}- (86) 
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(a) ferromagnetc core vortex 

Hp ~ H C2 



(0,h) 



(0,-W 



(b) polar core vortex 



Hp 



Hc 2 



(0,1 2 ) 

M 



(o,-i 2 ) 



► 



(0,1 2 ) 



(0,-l 2 ) 



(c (1) )* 
► 



H F 



(i F )* 



(<U 2 ) 



(0,-l 2 ) 



(0,1 2 ) 



(0,-l 2 ) 



H F 



(i F ) 



F'* 



H C2 



(Co H )* 



(0,1 2 ) 
— • — 



(0,-l 2 ) 



(0,1 2 ) 



(0,-l 2 ) 



Hp 



(ip)# 



(0,1 2 ) 



FIG. 5: (Color online) Deformation process for the case of (a) a FM-core vortex with n/ ■ — and (b) a polar-core vortex with 
n/ = 1 in the polar state, (a) The vortex, which is topologically stable in the polar state, shrinks to a point in the FM state. In 
the intermediate region, we run though the C2 state, (b) The vortex shrinks to a point in the polar state. In the intermediate 
region, we pass through the C2 state and the FM state. The FM state has a nontrivial winding. 



We call Eq. (|84|) a polar-core vortex, which was observed 
experimentally [||. 

Next, by using the map (c% )* (n p ^ 0), we can obtain 
different types of vortex-core structures. We obtain the 
following: 

(i F ).o(c< n ' ) ).o(ip).o(c(;)).([(o,-i 2 )]) 

= (*f). ° (4" p) )* ° (ipUKKnp-e-^™*)}) 
= (i F )*°(tf p) U[(™ p ,-e- i ^™*))) 
= (iF).([(7m p ,e- i ^)]) 
= [(7rn p ,e-^^)] 



i.e. 



[(0,1 2 )] eTTo^p), 



(87) 



where the last equality in Eq. (|57|) is obtained by the 
homotopy equivalence relations [see Eq. 



(0,l a )~(n p 7r,e- < * irn ')6(fl-p)o. 



The deformation process of Eq. (|87|l is illustrated in Fig.fJ] 
(b). Hence, Eq. (|87p implies that there is a vortex whose 
core structure is given by 



^hiroiHcJ^MHv), (89) 



(4" p) )* o (i P ). o (p£).([0,-l 2 ]) = [{im p ,e- iS ?™*% 

(90a) 

(i P ), o («£*>), o (» P ), o (c(;)),([(0, -l a )]) = [(0, 1 2 )]. 

(90b) 



We describe the order of strata as 



{(Vf> -+ (Vc 2 ) (</>p) (Vc 2 ) (V>f)}. (91) 



Since the core state is the FM state, we call Eq. (|89j) 
a FM-core vortex, which is considered as a new vortex- 
core structure. From Eqs. and (f59")l . we can achieve 

the vortex-core state when the map (cq ')* or (cg™ p ' ) )* 
appears in the sequence. Thus, if these maps do not 
appear, we can construct a long sequence of vortex-core 
structure such as 



}. (92) 



where this sequence ([BT)|) satisfies Eqs. (|53al) and (I53bl) . 



However, we do not discuss such structures further in 
this paper because the extension to this case is straight- 
forward. 

As a result, the core structure depends on the first 
choice of (dnj)*, and we can interpret n p as a quan- 
tum number of the core. In practice, n p is equivalent to 
the winding number along a ring of the polar state in the 
core. The FM-core vortex is accompanied by a surround- 
ing polar state with a nontrivial winding number n p ^= 
in the core. We call such a vortex, whose core structure 
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has a nontrivial winding number, an excited core state. 
In Sec.|VJ we show that a FM-core vortex is energetically 
stable under the quadratic Zeeman effect. 



2. Polar state 

For the case of c\ > 0, the boundary state becomes the 
polar state, 



^BS=Vp = t/^(l,0,l) J 



(93) 



Since itq(Hp) = Z, it is possible to create infinite 
types of vortices. The vortex is characterized by 
[(imp, e~ i ^ 7r ™ p )] S iro(Hp) (n p e Z). In what follows, 
we consider a vortex with the winding number n p . First, 

by using (cq" p ')», wc obtain 

= ( lF )M™ p ,e-^™?)}) 

= [(7rn p ,e-^^)] 

= [(0,1 2 )]. (94) 

The last equality in Eq. (|M|) is given by the homotopy 
equivalence relation given in Eq. ([55)1 . The deforma- 
tion process Eq. (|94|) is shown in Fig. [5] (a). Therefore, 
Eq. (|M|) means that there is a vortex whose core struc- 
ture is given by 

MHp) Tr (H c . 2 ) MHf), (95) 



where the sequence satisfies Eqs. (|53a[) and (|53b[) : 

(4 rlp) ),([(7rn p ,e-^^)]) = [(™ p , e"^™*)], (96a) 

(i F ). o (4"' ) ),(K ) e- 1 *^)]) = [(0, l a )], v n p e Z. 

(96b) 

The order of strata is given by 



(97) 



Accordingly, this type of vortex-core structure has an 
intermediate state involving the C2 state. By using 
(c^™ 1 ^)*, we obtain 



(i P ). o (4")), o ( iF ). o (c^>Ul(nn p ,e-^™r)}) 

= (i P )*o^-\o(i F U[(nn p -e- i ^)]) 

= (*p). o (^MK^-e"^™*)]) 
= (ip),([(0,-l)]) 

= [(0,-1)] 

= [(0,1)], (98) 

where the final equality is given by the homotopy equiv- 
alence relation in Eq. (|83|) . The deformation process in 



Eq. (|98|) is illustrated in Fig.[5](b). Hence, there exists a 
vortex-core structure given by 

7t (#p) ^hMH C2 ) ^ MHf) 



(ip) 



> MHc 2 ) MHp), (99) 

where sequence (|99| satisfies Eqs. (j53aj) and (|53h>() : 

° (<f). o (4 nf,) )*([(7rn p ,e-^^)]) = [(0,-l 2 )], 

(100a) 



(ip)* ° (4 _) )* o (z F ), o ( C i"^)*([(7rn p ,e-^^)]) 



(«p)i 



The order of strata is given by 



[(o,i 2 )]. 

(100b) 



(</>p)}. (101) 



Similarly to the FM state, we can prove the exten- 
sion of Eq. (|99|) to the long sequence vortex-core struc- 
ture, but we only state the result without discussing 
such a core structure: there are two types of vortex- 
core structures, which depend on which of (cg lp ' 1 )* and 

(cf 1 ^)* operates first. Similarly to the FM state, the 
label rif = {0,1} is interpreted as a quantum number 
characterizing the core, and it is equivalent to the wind- 
ing number along a ring of the FM state in the core. 
Therefore, there exists an excited state of the core in the 
polar state. 



V. ENERGETIC STABILITY OF VORTEX 
CORE STATES 

In this section, we discuss the energetic stability of 
a vortex-core state. We show that a vortex whose core 
structure has a nontrivial winding number can be realized 
experimentally. Wc consider a system with a negative 
quadratic Zeeman energy. The mean-field energy with 
the negative quadratic Zeeman term is given by 



E 



dr (e k i: 



Cint 



-zeeman 



(102) 



where ekin and ej n t are given by Eqs. (|9ap and (|9b|) . and 
the quadratic Zeeman energy density e zccman is given by 



-zeeman 



|«z|(h/>i| 2 + hM 2 )- 



(103) 



Here, q is the coefficient of the quadratic Zeeman term. 
The negative quadratic Zeeman is experimentally real- 
ized by the technique of microwave dressing [33j . 

To investigate the energetic stability of the vortex 
state, we solve the Gross-Pitaevskii equation and seek a 
stationary vortex state under a fixed boundary condition. 
We calculate the two-dimensional Gross-Pitaevskii equa- 
tion for spin-1 BECs in a uniform system with the Clank- 
Nicolson method. The time-dependent Gross-Pitaevskii 
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equations in the presence of a quadratic Zccman effect 



(104a) 



are given by 




dihi ( 

lh ~dt = \ 


h 2 

~2M 




h 2 

2M 


lh at = { 


h 2 

2M 



where (f z ), (/+), and (/_) are given by 



(/*) : = ^*n{fz)mni>n, 



(104b) 



(104c) 



(105a) 



</±> 



- y 



r m (f x ± »7v w», (105b) 



respectively. We seek stable solutions by an imaginary 
time propagation with the fixed boundary condition. We 
assume cq 3> \c\\. Therefore, the maximally enlarged 
OPM is S 5 , which is consistent with our classification 
method discussed above. 

For the FM state (ci < 0), we choose the boundary 
condition as 



t/> F = V^,o,o) 5 



(106) 



where <p is the azimuthal angle in the two-dimensional 
space. The numerical results are shown in Fig |5J where 
we show stable solutions for q = and q < q c , respec- 
tively. Here, q c is a transition point from the polar-core 
vortex to the FM-core vortex. At q = 0, the vortex core 
is filled with the polar state as shown in Fig. [5] (b). The 
vortex-core state is given by ipp = (0, y/n, 0) T . Since 
the intermediate region is covered by the C2 state, this 
is the polar-core vortex whose vortex-core structure cor- 
responds to Eq. (j54"| and the order of strata is given by 
Eq. (|86| . The polar-core vortex has been realized exper- 
imentally @, and our classification method is consistent 
with both the experimental observation and the numeri- 
cal simulation. 

On the other hand, at q < q c , the vortex core is filled 
with the FM state. As shown in Fig.[6](c), we have a ring 
of the polar state in the intermediate region. Because the 
component ipo is negligible for q < q c due to -Ezocman oc 
IgllV'ol 2 - Hence, the core state becomes = (0, 0, ^/n) T 
and the polar state on the ring is given by 



(107) 



Here, Eq. (|107p represents the half-quantum vortex. 

Therefore, the core structure corresponds to Eq. 

and the corresponding order of strata is given by Eq. (|91[) 



(a) 
FM state 



C2 state polar state 



• • t 

(b) polar-core vortex 



(c) FM-core vortex 



Half-quantiezed vortex 




(d) magnetization 
polar-core vortex 



FM-core vortex 




FIG. 6: (Color) (a) Plots of \ip{fi,(j>)\ 
the C2 state, and the polar state, 



for the FM state, 
where %l)(6,cj>) = 



X)m=-i Yim(0, (j))ipm- Here, Y\ m is a rank-1 spherical- 
harmonic function and colors represent the phase of if)(8,(f>). 
(b) and (c) Vortex-core structures of the polar-core vortex 
(kl/|ci| = 0) and the FM-core vortex (|g|/|ci| = 0.4) by the 
spherical-harmonic representation. We assume co/|ci| = 100. 
The polar-core vortex is consistent with the core structure in 
Eq. (|84[) . whereas the FM-core vortex corresponds to the core 
structure in Eq. (|89|) . In the intermediate region, there exists 
a half-quantum vortex along a ring of the polar state, (d) 
Amplitude of the magnetization \F\ in the two-dimensional 
system. The vertical and horizontal axes show the scaled di- 
mensionless y and x axes, respectively. The red and blue re- 
gions represent the FM state and the polar state, respectively. 
The other colors show the C2 state. 
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with rip = 1. Since the vortex core has nontrivial wind- 
ing, it is an excited state of the core. 

Next, in the polar state (ci > 0), we choose a boundary 
condition as the half-quantum vortex in Eq. (|107|) . For 
q < 0, the vortex-core is always filled with the FM state 
i/>f = (0, 0, \/n) T . The order parameter in the intermedi- 
ate region is always the C 2 state. Hence, the vortex-core 
structure corresponds to Eq. (|9"5")l and the order of strata 
is given by Eq. (|97|) . This vortex-core structure has been 
discussed numerically in Refs. 0] and [T^]. The excited 
core state is unstable at least in the presence of a negative 
quadratic Zeeman effect. 



VI. SUMMARY AND DISCUSSION 

In this paper, we have developed a method to classify 
vortex-core structures. Our classification generalizes the 
previous works 0, HH so as to identify the state that 
fills the vortex core when the vortex is nonsingular. Our 
classification method is formulated as follows: first, we 
define the maximally enlarged OPM, which is determined 
by a constant particle density, and decompose this max- 
imally enlarged OPM into strata. By defining the map 
between the strata, we have characterized the vortex-core 
structure. 

We have applied this method to the spin-1 BEC, and 
as a result, wc have found that the vortex-core structure 
is made up of several states which are aligned in a con- 
centric fashion and, as a whole, becomes an excited state 
of the core. Each state of the concentric ring has a local 
winding number which constitutes a quantum number 
characterizing the vortex-core structure. 

For the case of the FM state, the vortex-core structure 
is formed from the concentric pattern of the FM state, 
the C 2 state, and the polar state. Under the minimal 
concentric pattern, the vortex-core structure is classified 
by an integer, which represents a winding number de- 
fined along a ring of the polar state in the core. For the 
case of the polar state, the vortex-core structure is also 
composed of the FM state, the C 2 state, and the polar 
state. There exists one nontrivial vortex-core excitation. 
We have numerically examined the energetic stability of 
the vortex-core structure that has a nontrivial winding 
number and found that it can be stabilized in the FM 
state under the negative quadratic Zeeman effect. 

In this paper, we have focused on the spin-1 BEC; 
however, our method can also be applied to the case of 
higher spins. These applications will be discussed else- 
where. We comment here on the role of linear Zeeman 
shift. Since the cold-atom systems always conserve the 
total spin, the linear Zeeman shift merely induces the 
Larmor precession, which can be eliminated by moving 
onto the rotating frame of reference in spin space. How- 
ever, if we consider the case in which the total spin is 
not conserved, we cannot ignore the linear Zeeman shift. 
In this case, the direction of the magnetization becomes 
parallel to the magnetic field, and therefore the order pa- 



rameter is described by a single component. Hence, the 
maximally enlarged OPM is S 1 . Since tti(S 1 ) = Z, the 
vortex core is singular. 

Finally, we discuss yet another application of our clas- 
sification method to the stability of a vortex under the 
transition between two ordered phases. We can easily 
modify our scheme from the classification of the vortex 
core to the stability of a vortex under the transition by 
replacing distance r with time t. If wc apply our method 
to the stability of a vortex under the transition, there ex- 
ists another deformation process because the initial state 
need not have the winding number. For example, when 
the initial state is the FM state and the final state is the 
polar state, we run though the C 2 state at an intermedi- 
ate stage. It is possible to change to the C 2 state with 
a nontrivial winding number n p even when the FM state 
is uniform. 

Note added in proof Recently, a related work by 
J. Lovegrove, et al. [34| appeared which discuss as ener- 
getically stable vortex state in a spin-1 BEC. 



ACKNOWLEDGMENTS 

The authors would like to thank D. A. Takahashi and 
M. Kobayashi for useful discussions. This work was sup- 
ported by Grants-in-Aid for Scientific Research (Kakenhi 
Nos. 22340114, 22103005, 22740265, and 23740198), a 
Grants-in-Aid for Scientific Research on Innovative Areas 
"Topological Quantum Phenomena" (Nos. 22103005 and 
23103515), a Global COE Program "The Physical Sci- 
ence Frontier" , and the Photon Frontier Netowork Pro- 
gram of MEXT of Japan. S.K. acknowledges support 
from JSPS (Grant No.228338). Y.K. acknowledges the 
financial support from Inoue foundation for science. 



Appendix A: CALCULATION OF EQ. (^5j) 

A spin-1 



Let us describe the calculation of Eq. (|55|) . 
parameter is given by 

t/> = (V>i,Vo,V>-i) T e C 3 , 

where ip m (m = 1, 0, —1) are described by 

tpi = x\ + ix 2 , 
4>o = x 3 + 2X4 , 
ip-i = x 5 + ixq. 

Here, x { 6l(i = l,2, 3, 4, 5, 6). To derive Eq. 
operate e l ^ v ^e l ^ a on ip: 



(Al) 

(A2a) 
(A2b) 
(A2c) 

(jST?)) . wc 



e"*cos 2 § ^ 

sin 8 

^7T 



Via • 2 8 sin a 
e sin f - ~jf 



e~ ia sin 2 4 



,-iq sin p 
' la COS 2 f 





h ix 2 \ 


)(«: 


- 1X4 




h iXQ J 



(A3) 
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We determine a and (3 so that the m = component 14 ^ as 
vanishes. Hence, we obtain a and f3 for i 3 ^ and 



a = cos 



(3 = cot" 1 



xa{x 2 + xe) + x 3 (xi + x 5 ) 



yj {Xi{x b - Xi) - X 3 (x 6 - X 2 )} 2 + {Xi{x 2 + Xfjj + X 3 (xi + X 6 )} 2 

I 'y>2 m2 I /y>2 2 

1 x 5 ^1 ' x 6 x 2 

V2 ^ {x 4 (x 5 - sci) - x 3 (x 6 - x 2 )} 2 + {xi(x 2 + x 6 ) + x 3 (xi + x 6 )}< 

I 



(A4a) 
(A4b) 



Therefore, the remaining components are described by 
x'i, x'2, x' 5 , and x'§ 



X — \~ ix^ 



x' 5 + ix' 6 



(A5) 



Next, we operate e^^e**" 7 on Eq. (|A5[) and we give <j> 
and 7 in order for the remaining components to be real 
as 



7 = 



2 

h + O-i 



where 9\ and 8-1 are defined by 



h = tan _± -f 

-1 / ^6 



-1 =tan^ -± 



(A6a) 
(A6b) 

(A7a) 
(A7b) 



Giving a, /3, 7, and as Eqs. (|A4a|) . (|A4bl) . 
(|A6a[) and (|A6b() . the remaining components become 

y/x[ 2 + x' 2 2 and ^Ac' 5 2 + a/ 6 2 . Using Eqs. (HUJ) and dHJ, 
the order parameter is given by using n and |f| as 

( 



Thus, we obtain Eq. (|55|) . 



(A8) 



map pi is defined by the map from G to G such that 

Pi : G -> G, 
(x, e"^ 8 ) ^ (e", e"^ 9 ) v x el, v fie [0, 4tt), (Bl) 

where v = x,y, z, cr„ is a Pauli matrix, and pi is a sur- 
jective and homomorphic map. Hence, since e l27m = 1 
and e l f" 2 '* = I3, we can obtain a kernel of pi as 

kerpi = {(27m, 1 2 ), (27m, -l 2 )\n g Z}, (B2) 

Second, P2 is defined by the map from G to G/H q such 
that 

p 2 : G -> G/JT„ 
( e #, e -'/"«) ^ [( e <* e"^ Q )] ? v 0, a e [0, 2tt), (B3) 

where [• • • ] q represents an equivalence class defined by 

gg'- 1 eH q & ~ g' for any g,g' g G. (B4) 

Here, p 2 is also a surjective and homomorphic map. By 
definition (|B2j) . we obtain the kernel of p 2 by 



kerp 2 = ff,. (B5) 

By using pi and p 2 , we make up the map from G to G / H q 

P2 o Pl :G^G/H q . (B6) 

Thus, H q is rewritten as 

H q = kcrp 2 op 1 . (B7) 

To calculate we define the inclusion map from G to 
G as follows: 



Appendix B: CALCULATION OF H F , H F , AND Hc 2 

In this section, we give the calculation of the lifting 
isotropy group for each G-orbit. Each isotropy group 
is given by Eqs. (fT5]l. (fT5J). and ([55]). The universal 
covering space of G is given by Eq. (p2"j) . To calculate 
i/g, we define two projection maps pi and p 2 . First, the 



t (0lla) :(e^,e~^)^(0,e" 1 *''), 

^(27r,l 2 ) 
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Then, we can give the lifted isotropy group H q by using We first calculate Hp. Since H-p is given by Eq. (IT8l) . H-p 
l s (s € hevpi) as is obtained as 

H q =p^\H q ) (B8) 

= (J L(2-Kn,l 2 )(h) U (J t(27rn,-l 2 )(^)i ^ € 

(B9) 



H F = (J t(2 ™,i 2 )((e ie ,e-^ e ))U |J t (27rri ,_ l2) ((e*V-^ )), £ [0,2tt), 



= • • • U (0 - 2tt, e"^0) U (0, e~^) U (0 + 2tt, e^ 1 ^) U • • • U (0 + 2im, e"* = y ) U • • • , 
• • • U (0 - 2tt, -e _i ^0) U (6», -e~ ?; * e ) U (0 + 2tt, -e" 1 * 6 ) U • • • U (0 + 2nn, -e~' l ^ e ) U • • • , 
= • • ■ U (0 - 2tt, e-iTfP-a*)) u (0, e^ 9 ) U (0 + 2tt, e -^(e+2^)) y . . . j 
■ • • U {9 - 2tt, -e-**^- 2 *)) U (0, -e-^ e ) U (0 + 2tt, -e"^^ 2 *)) U • • • , 

= (x,e ts **)U (x,-e iS * x ), (BIO) 

I 

where xgl and e ±l ^- 2lT = — 1 2 . Next, we calculate Hp. where n € Z and e [0,47r). Finally, we calculate Hc 2 - 
Since -ff P is given by Eq. ([19]), ff P is given by Since Hq 2 is given by Eq. ([59]), H C2 is given by 

ffp = U ^n,±x 2) {{he-^)) (Bll) ^ = y t(27rn±l2)((1)l3))u (J t(2 ^ ±l2) (( e -,e-/-)), 



U U Wi.)^^.^- 4 '-)). = K±l2 )u(( 2li + lK ± e- 

(27rn, e'^ 8 ) U ((2n + 1)tt, e^V 4 *") , 



(B13) 



where n e Z. Therefore, Hc 2 is composed of discrete 
(B12) 

elements. 
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